Nonlinear effects of phonon fluctuations on transport through nanoscale junctions 
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We analyze the effect of electron-phonon coupling on the full counting statistics of a molecular 
junction beyond the lowest order perturbation theory. Our approach allows to take into account an- 
alytically the feedback between the non-equilibrium phonon and electronic distributions in the quan- 
tum regime. We show that for junctions with high transmission and relatively weak electron-phonon 
coupling this feedback gives rise to increasingly higher nonlinearities in the voltage dependence of 
the cumulants of the transmitted charges distribution. 
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Single molecule junctions and atomic chains suspended 
between metallic electrodes constitute a fascinating play- 
ground to explore the interplay between electronic and 
vibronic degrees of freedom, sec e.g. Rcfs. [lH3. The 
interest is now not only restricted to the understand- 
ing of the mean current-voltage characteristics but has 
been extended to noise properties^ and, more generally, 
to the full counting statistics (FCS) of the transmitted 
chargesi^ An intense theoretical activity has been fo- 
cussed on the analysis of the simplest model consisting 
in a resonant level coupled to a single phonon mode in the 
quantum coherent rcgimci^r— So far, however, several as- 
pects of this problem remain to be clarified. A serious 
limitation of existing transport theories is that they do 
not take into account the influence of the non-equilibrium 
phonon fluctuations in the statistics of the transmitted 
electrons, namely the feedback of the phonon dynamics 
on the current-noise properties'^ This limitation is as- 
sociated with the breakdown of perturbation theory be- 
yond the lowest order in electron-phonon (c-ph) coupling 
as reported in Refs. [31 . These works, which demon- 
strate the necessity of including non-perturbative effects 
in the analysis, are, however, limited to equilibrium prop- 
erties of bulk materials or to individual molecules in the 
sequential tunneling regime. 

In this work we demonstrate, by a partial resummation 
of the perturbative expansion, the great impact of the 
feedback of the phonon dynamics on the quantum trans- 
port properties through nanoscale junctions with high 
transmission and relatively weak e-ph coupling. The ac- 
tual signatures of phonon fluctuations result from the 
interplay of several energy scales, i.e. the tunneling 
rate T of electrons, the e-ph coupling A, the phonon 
frequency ujq and the relaxation rate oc 77 of the local 
phonon mode due to the coupling with bulk plionons. 
Depending on the specific sample considered and the effi- 
ciency of the relaxation mechanism for the phonon popu- 
lation, one might obtain a regime characterized by a ther- 
mal phonon population when rj 3> A^wo/r^ (equilibrated 
phonons) or a regime where a strong non-equilibrium 
population is generated when rj <^ A^wq/T^ (unequili- 
brated phonons). We analyze the crossover between the 
two regimes and demonstrate that in the regime of un- 
equilibratcd phonons the electronic current-noise shows 



a strong nonlinear behavior as a function of the applied 
voltage for V > luq. We attribute these nonlinearities 
to a dynamical feedback of electronic quantum fluctua- 
tions that strongly rcnormalizc the parameters describing 
the local phonon mode. This mechanism affects also the 
higher-order cumulants C„ of the current operator which 
exhibit growing nonlinearities with the order of the con- 
sidered cumulant. In particular we find the scaling 



SCn+l/SCn ^ V/UJO 



(1) 



for the A-dependent part (5C„ of the cumulants which is 
valid for r V ^ ujQ ^ V- Despite the experimen- 
tally weak e-ph coupling these strong nonlinearities be- 
come sizable in the regime corresponding to recent exper- 
iments on transport through small molecules using me- 
chanically controllable break j unctions ii Therefore they 
are expected to be essential to correctly capture the noise 
properties above the inelastic threshold. 

For our study, we use the model sketched in Fig. [TJa) 
corresponding to a Hamiltonian 



H = AdU + cjofl^a + Hl + Hb + Vt + Ve- 



-ph • 



(2) 



Here d is the annihilation operator for the electron on 
the dot level at energy A and a is the annihilation op- 
erator for a phonon of energy luq (we use e = ^ = 1). 
The left and right electrodes, described by H]^/fi, are 
modeled as noninteracting fermionic continua with re- 
spective field operators ^'l/aC^;) and the corresponding 
electronic density of states po is assumed to be energy 
independent. The electrodes arc held at different chem- 
ical potentials giving rise to a voltage drop V > 0. 
Electrons are allowed to tunnel between the electrodes 
and the dot, described by the local tunnel Hamiltonian 
Vt = d1'(7L^'L(0)+7i7*i7(0))+h.c. with left(right) hop- 
ping amplitude ^l{r)- Finally, the e-ph interaction is 
given in terms of the displacement operator $ = a + 
by Ve~ph = X^d'^d. In the present work we concentrate 
on the regime Tlh = t^PojI ^ y,a;o,A,A which is 
relevant for experiments involving small molecules^ and 
chains of atoms 1^ The analysis will mainly address the 
case of perfect transmission (F^ = Tn = F/2 and A ~ 0) 
where the feedback mechanism can be explained in terms 
of simple analytical expressions. 
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FIG. 1: (a) Sketch of the model system, (b) Exchange(Fock) 
and (c) Hartree diagram, both appearing in the linked cluster 
expansion of the CGF at order A^. Bottom: fourth order 
diagrams, i.e. (d) double- bubble, (e) double-exchange, and 
(f) rainbow diagram. 



The FCS of charge transport through the dot is derived 
by calculating the cumulant generating function (CGF) 
xiO = XoiO + SxiO of the corresponding probability 
distribution P{Q) to transfer the charge Q during the 
measuring time <o- The CGF of the non- interacting sys- 
tem xoiO is given by the Levitov-Lesovik result.— The 
correction due to finite c-ph coupling is given by 



MO ^-\^{Tc CXp -i / dtVe-ph{t)dt 



(3) 



Here the expectation value is taken with respect to the 
non-interacting Hamiltonian Hq = H — Ve-ph and inte- 
gration runs over the Keldysh contour C with time or- 
dering Tc along C. The counting field ^{t) is included 
in the tunnel Hamiltonian Vr through the substitution 
7i = 7^6^''^'*^/^, taking the value +(— )C on the for- 
ward(backward) branch of C. The CGF allows to access 
the cumulants C„ of P{Q) by successive derivation with 
respect to the counting field, C„ = f|^|g=o- 

The linked cluster expansion of 5x in the e-ph coupling 
A can be evaluated in terms of the noninteracting electron 
Green function (GF) Ga{t,t') = -~i{Tcd{t)d) {t'))^ and 
phonon propagator Do = -i{Tc<^{t)^{t')) ^. The 
Keldysh components (a,/3 = ±) of the latter read 



1 



1 



aijjQ + if] u) + l3ujo — 11] ' 



(4) 



which includes the finite phenomenological parameter 
77 describing the coupling with bulk phonons. The 
first term contributing to the expansion of Sx is c>c 
and involves the calculation of a Hartree-type and an 
exchange(Fock)-type diagram, c.f. Fig. [Ub)-(c). The 
latter term was recently discussed in detail^i^ and was 
shown to be responsible for the jump in the derivative 
with respect to voltage of the cumulants near the inelas- 
tic threshold at V ^ ljq. 



On the other hand, the Hartree diagram leads to a 
shift of the dot level position A and this renormalization 
is almost independent of voltage. This statement is found 
to hold also for the Hartree-type diagrams of higher order 
in A within the parameter range discussed in this article. 

The next order term A'*(Sx4(^) involves the evaluation 
of three diagrams, which are shown in the lower panel 
of Fig. [TJ The iromenclature double bubble (DB), double 
exchange (DX) and rainbow (RB) refers to the Feynman 
graphs of the corresponding self-energies. In the limit of 
r ^ V, cjQ and at zero temperature, we can give analyt- 
ical expressions for the DX and RB diagrams and find 
that they are of order A^'/F^. They can thus be safely 
neglected in comparison to the dominant DB diagram 
Sx4^^\0 oc A*/(r^ry^) which diverges in the limit of 
vanishing external damping 77. 

With the corresponding diagrams of higher order in 
A all being found to be divergent, we suggest a rcsum- 
mation of this subset of diagrams. The Dyson equation 
D = Dq + DoTjphD for this RPA-typc dressing of the 
phonon GF involves the ^- and T^-dependent phonon self- 
energy 
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n++(cj) -H+-(a;) 
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where n(w) is the charge polarization loop, 

n"^(o;) = J ^Gf{u^ + e)G^,'^{s). (6) 

Note that in the parameter range discussed in the follow- 
ing, i.e. A^/wqF < 0.25, we found the RPA approxima- 
tion to be charge conserving up to the chosen numerical 
accuracy. The main result of the rcsummation is a shift 
and a broadening of the phonon mode by adding a finite 
lifetime due to the coupling to the electrons. A crucial 
aspect is to realize that this mode renormalization is 
dependent, i.e. it is sensitive to the counting-field of the 
tunneling electrons. This has important consequences as 
discussed below. In the limit T ^ V,u} we can expand 



L cr=± 



. " 2\u\ + {V-\i,\)e{V-\^\) 
" ^""'-JnT^ 2nr^ 

and thereby find four poles fi^ of D"^{uj). For V > ljo, 
the two poles having a positive real part read 



A^ A^ 

ni/2 = {OJQ ±177) — - 



ttF 27r2F2a;o 



(9) 



and due to symmetry ^3/4 = — ^1/2. 

This phonon renormalization goes beyond the widely 
used approximation of the non-equilibrium phonon pop- 
ulation based on a rate equation (RE) descriptioni^i^ 
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FIG. 2: (a) Mean number of phonons (n) a.t V — SuJo and 
cjo/r = 0.01 as a function of 77 shown for various values of A. 
For comparison, the dotted curves show the results of the rate- 
equation approach, Eq. (|10|l . for the same parameters, (b) 
Voltage dependence of (n) for 77/r = 10~^ and too/T = 0.01. 



FIG. 3: (a) Differential conductance and (b) derivative of the 
current noise with respect to voltage as function of voltage 
for various values of A. Parameters used are loq/T = 0.01 
and 77/r = 10~^. The dotted (dashed) lines show the results 
of the rate equation approach for thermally equilibratedfun- 
equilibrated) phonons for comparison, taken from Ref. 



Within this RE approximation, the mean number of 
phonons (n) at zero temperature is given hy^ 



(n) 



RE 



UJQ 



1 



\V\-uo), (10) 



which quahtatively describes the crossover between the 
equihbratcd and unequilibrated phonon regimes (for 77 ^ 
A^wo/r^ and 77 < A^wo/r^, respectively). In Fig. [IJa) 
we compare this prediction (dotted hnes) with the re- 
sults for (n) obtained within the renormalized pertur- 
bation theory (solid lines) for luq = O.Oir and different 
A values. Therefore we compute (77,) by integrating the 
non-diagonal component of the dressed phonon propa- 
gator, (n) — — i + ^/ duj D~^~ {uj)\^=o- While there is 
a good agreement of the RE approach with our calcu- 
lations when A < wq, discrepancies arise when A > wq, 
especially in the regime of unequilibrated phonons. The 
RE approach then underestimates the phonon popula- 
tion at a given voltage due to the neglect of the mode 
renormalization effects. 

We can now turn to the calculation of transport prop- 
erties. We evaluate SxiO using the dressed phonon GF, 

aJ3=± '' " 

(11) 

and from this deduce the current-voltage characteristics 
I{V) as well as the current noise power (second cumu- 
lant C2). Figure [3] shows representative plots for the 
A-dependent contributions 5G and 5S' to the differential 
conductance G = dl/dV and derivative of the current 
noise with respect to voltage S' = dC^/dV , respectively, 
in the crossover regime 77 ^ A^wo/r^. The obtained 



curves change dramatically with the e-ph coupling. The 
conductance (Fig. [5^) shows a negative jump at ^ = wq 
discussed in detail in recent works^"— which is related 
to enhanced inelastic backscattering processes. Note 
that the rounded spikes at the threshold are reminiscent 
of the log-singularities obtained without the resumma- 
tion which are pronounced only for the largest A val- 
ues shownii^ For voltages above the inelastic threshold, 
5G{V) is almost linear and exhibits a crossover from a flat 
voltage-independent regime at low A/F — 0.001 (equili- 
brated phonons) to a negative slope at larger A/F = 0.05 
(unequilibrated phonons). The corresponding limiting 
cases are plotted (dotted and dashed hnes, respectively) 
using the RE approach of Ref. and the shape of the 
5G{V) curves is well reproduced by this approximation. 
The linear voltage dependence is directly related to the 
linear population of the phonon mode for V > ojq, c.f. 
Fig. [2Ib) , and the negative slope is due to activated ab- 
sorption processes resulting from the increasing popula- 
tion of the local vibrational mode. 

The derivative of the current noise with respect to volt- 
age (Fig. [SJd) exhibits a positive jump ai V = ujq as. & 
result of the same inelastic backscattering mechanism re- 
sponsible for the jump in 5G{V). However, for V > ujq 
the behavior of the curves is qualitatively different. In the 
regime of equilibrated phonons, 5S' {V) is almost voltage- 
independent and well reproduced by the RE approach. 
For unequilibrated phonons however, the current-noise 
curves develop strong nonlinearities which are absent in 
the RE approach. We stress that already for quite small 
A/F = 0.005 (corresponding to an almost thermally equi- 
librated regime, c.f. Fig.[2]a,) the amplitude of the nonlin- 
earities is large enough to match the (maximal) results 
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of the RE approach for unequihbrated phonons (dashed 
hne in Fig.[3]D). This quahtative different behavior can be 
understood as follows. Within the RE approach and in 
the limit F 3> u!q,V,X in which we are interested, the 
correction to the FCS due to e-ph coupling is purely 
PoissoniaUji^i^, i.e. SCn.RE oc SCi.re- On the other 
hand, the correction to the FCS in the rcnormalizcd the- 
ory departs from being Poissonian due to the dynamical 
renormalization of the phonon mode. In the particular 
regime of unequihbrated phonons (77 <C A^w/F^) and for 
F ^ F > cjq, the dressed phonon propagator can be 
approximated by 



D 



n 



with f = {lu'^ — ijS^)/2uj^ and the poles Vlj given in Eq. 
dni). The integral of Eq. pT|) is then approximated by 



2-K 



(13) 



where the prefactor can be expressed as 
4a;g zttF^ e'^ ' 



itan(^/2) 



V 



too 



(14) 



with aj = — sign(Im(J7j)). This factor therefore is identi- 
fied as the source of the increasing nonlinear voltage de- 
pendence when taking successive derivatives of Sx with 
respect to ^ at ^ = 0. Finally, the most nonlinear contri- 
bution to the SCn in this regime reads 



n+1 



(15) 



which demonstrates the scaling law anticipated in the 
introductionJi 

In conclusion, we have analyzed the effect of c-ph in- 
teraction on the FCS of a molecular junction beyond 
the lowest order perturbation theory. A RPA-type re- 
summation of diverging diagrams results in a renor- 
malization of the phonon parameters which depends 
on the electron counting- field, thereby giving rise to 
an increasingly higher nonlinear voltage dependence in 
the cumulants of the transmitted charges distribution. 
These strong nonlinear effects are expected to be accessi- 
ble in low-temperature transport experiments performed 
on highly conductive molecular junctionsi and atomic 
chains?^ where typically the energy of the strongly cou- 
pled vibrational mode and the e-ph coupling strength are 
of order wq/F « 0.01 and A/F ss 0.1 — 0.01, respectively. 
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